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Abstract. We experimentally study the transient motion of a colloidal particle 
actively dragged by an optical trap through different viscoelastic fluids (wormlike 
micelles, polymer solutions, and entangled A-phage DNA). We observe that, after 
sudden removal of the moving trap, the particle recoils due to the recovery of 
the deformed fluid microstructure. We find that the transient dynamics of the 
particle proceeds via a double exponential relaxation, whose relaxation times remain 
independent of the initial particle velocity whereas their amplitudes strongly depend 
on it. While the fastest relaxation mirrors the viscous damping of the particle by the 
solvent, the slow relaxation results from the recovery of the strained viscoelastic matrix. 
We show that this transient information, which has no counterpart in Newtonian 
fluids, can be exploited to investigate linear and nonlinear rheological properties of 
the embedding fluid, thus providing a novel method to perform transient rheology at 
the micron-scale. 
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1. Introduction 

Viscoelasticity is ubiquitous in many materials ranging from biological fluids, polymers, 
micellar systems, colloidal suspensions and more. In general, such materials are 
characterized by a strongly time-dependent mechanical response to stress or strain. 
For instance, in case of an oscillatory shear, they may exhibit either liquid or solid- 
like properties depending on the imposed frequency [1]. Such non-Newtonian behavior 
originates from the storage and dissipation of energy within their complex microstructure 
which gives rise to a hnite macroscopic stress-relaxation time [2] . Given the abundance of 
systems which involve the interplay between microscale viscoelastic flows and embedded 
micron-sized objects, as typically found in lab-on-a-chip devices 131 H i , porous 
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media EE], and in active matter [HiEiiin], the understanding of the equilibrium and 
non-equilibrium dynamics of colloidal particles in viscoelastic fluids is a noteworthy topic 
in soft matter. Microrheology has become in recent years a held which successfully 
addresses some of these problems. In particular, it provides alternative methods to 
bulk rheology to investigate how and deformation properties of microlitre samples of 
viscoelastic materials. For example, linear shear moduli can be inferred in passive 
microrheology from the thermal huctuations of suspended colloidal particles via a 
generalized Stokes-Einstein relation [m na [la [n]. On the other hand, in active 
microrheology small-amplitude oscillatory forces can be applied to the particle, which 
enables the dynamical measurement of the linear viscoelastic response of the surrounding 
huid [T^ [TB| [T71ITB] . More recent developments of active microrheology have aimed to 
probe nonlinear response of complex huids by perturbing their microstructure far away 
from equilibrium. This has been achieved by dragging or rotating the probe through 
the huid at sufficiently large constant force or torque [ISl [20112111221123] , or at constant 
velocity [2a ESI EH EZ], thus inducing e.g. thinning. 

Despite multiple potential applications of microrheology, so far most of the 
experiments have focused on the study of linear-response or steady-state quantities. 
More complex transient behavior observed macroscopically, e.g. nonlinear creep, 
stress relaxation and strain recovery, is far less well understood within the context 
of microrheology. Indeed, theoretical works suggest that the drag force acting on a 
colloidal probe driven through a non-Newtonian huid, strongly depends on the spatio- 
temporal deformation of the microstructure of the surrounding huid [291 Eol EH- 
Therefore, in contrast to Newtonian liquids, where a low Reynolds number guarantees 
simple Stokesian hows and linear-response drag forces, it is questionable whether this 
is applicable when a particle is dragged by a time-dependent driving force over large 
distances through a viscoelastic medium [32] • Then, it is not evident which kind of 
rheological properties of the huid can be extracted by simply tracking the position of a 
probe in response to a pulsed or Heaviside time-dependent driving force. It should be 
mentioned that such a time-dependent force represents the microrheological analogue 
to a macroscopic how startup/cessation experiment [2H], where the time evolution of 
the stress (strain) is measured upon suddenly imposing or relieving strain (stress), 
respectively. Although relevant for further developments of microrheology and for the 
general understanding of the motion of micron-sized particles in viscoelastic media, these 
questions have not been experimentally investigated so far. 

In this work, we study the local microstructural recovery of several viscoelastic 
huids by means of microrheology. For this purpose, we drag spherical particles at 
constant velocity through them using optical tweezers, thus inducing strain, and then 
we studdenly remove the trapping force. We observe that the particle recoils after 
removal of the trap until the complete recovery of the investigated huid. We hnd that 
the transient particle dynamics during the recovery proceeds via a double exponential 
relaxation, whose relaxation times remain independent of the initial particle velocity 
whereas their amplitudes strongly depend on it. These processes reveal how the elastic 
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Figure 1. (a) Experimental protocol to study the transient motion of a colloidal probe 
during the recovery of a viscoelastic fluid; (i) passive, (ii) active, and (ill) transient 
microrheology, (b) Trajectory of a r = 1.6 /rm particle suspended in a wormlike micellar 
solution (CPyCl/NaSal 7mM, T = 303 K) during the microstructural recovery after 
removal of the trap ik = 8pN/im“^, v = 3.20/rms“^). The blue arrow indicates the 
direction of the particle recoil. See text for explanation. 


energy stored by the fluid due to a large local strain is eventually dissipated. While 
the fastest relaxation mirrors the viscous damping of the particle by the solvent, the 
slow process results from the relaxation of the viscoelastic matrix. In particular, we 
show that this experimental method allows to probe linear and nonlinear transient 
rheological behavior of the viscoelastic fluid under study and to extract unambiguously 
its stress-relaxation time. In addition, our findings have important consequences on 
the understanding and interpretation of the motion of colloids actively driven through 
viscoelastic media. Thus, we provide clear experimental evidence that, unlike passive 
and small-amplitude oscillatory microrheology, the active motion of a colloidal particle 
over large distances does not only probe the bulk viscoelastic properties of the fluid but 
is also affected by the flow field and deformation imposed by the particle. 

2. Experimental description 

In our experiments, we prepare equimolar solutions of cetylpyridinium chloride (CPyCl) 
and sodium salicylate (NaSal) in deionized water at concentrations between 5 and 9 mM. 
We also studied samples of other viscoleastic fluids, namely, aqueous polyacrylamide 
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(PAAm) solutions (M^, = 18 x 10® at 0.05 % wt) and entangled A-phage DNA solutions 
(500 /ig/mL in 10 mM Tris-HCl and 1 niM EDTA). We add to each viscoelastic solution a 
small volume fraction of spherical silica particles, whose radii are chosen between r = 0.9 
and 2.2 /rm. The temperature T of each sample is kept hxed by a flow thermostat 
between 20 ± 0.1 and 40 ± 0.1°C. Then, a single particle is trapped by optical tweezers 
created by deflection of a Gaussian laser beam (A = 1070 nm) on a galvanostatically 
driven pair of mirrors and subsequent focusing by a microscope objective (100 x, NA 
= 1.3) into the sample. In order to avoid hydrodynamic interactions, the particle is 
trapped at least 20 fim away from the walls of the sample cell. The optical trap at 
position To = {xo,yo) exerts a restoring force on the particle at position r = {x,y), i.e. 
fo = —/c(r — To). Using video microscopy, we track the center of mass of the particle at 
120 frames s“^ with a spatial accuracy of 10 nm. We implemented the protocol sketched 
in Fig. [11(a), where the red solid line represents the time evolution of Xo, to study the 
recovery of the different viscoelastic fluids. 

(i) During time to < t < ti, the trap is kept at rest {xo = 0, |/o = 0), so that the particle 
is in thermal equilibrium with the surrounding fluid. The thermal fluctuations of 
{x, y) enable us to determine the stiffness k of the optical trap and to characterize 
the linear microrheological properties of the fluids by passive microrheology. 

(ii) Then, at time t = ti, we start to move the trap at constant velocity v along a linear 
path, {xo = v(t — ti), yo = 0), thus actively dragging the particle through the fluid 
along X and inducing strain. At position Xo = Xq, where Xq = 16 /rm r, the 
motion of the trap is reversed: v —)■ —v. We perform a sequence of these back-and- 
forth movements of the trap position until the distance travelled by the particle 
during a period Tq = 2Xofv remains constant. Depending on the value of the 
driving period Tq we can induce either linear or nonlinear microrheological response 
of the viscoelastic fluid |Z7|. During this step, we measure the effective viscosity of 
the fluid 7] = fd/{QTirv) by active microrheology from the balance between the drag 
force fd and |fo|. 

(hi) Finally, at time t = t 2 and position Xo = XqI2 the optical trap is suddenly turned 
off. During this hnal step, the particle motion at time t > ^2 is subjected to the 
relaxation of the previously strained microstructure of the fluid. The subsequent 
recovery results in a recoil of the particle along the x-direction opposite to that 
of the prior trap motion. This recoil is illustrated in Fig. [I](b), where we plot an 
example of a 2D trajectory of the particle during the recovery process of a wormlike 
micellar solution after removal of the moving trap. We observe that, in absence of 
external forces, the particle gradually performs Brownian motion once the fluid 
completely equilibrates. 

Note that, in a Newtonian fluid, the particle undergoes Brownian diffusive motion 
immediately after removal of the force and therefore such a recoil is not observed, as 
verihed in our experiments in H 2 O and glycerol (data not shown). Since the recoil takes 
place in the direction parallel to the driving force, in the following we only concentrate 
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Figure 2. (a) Time evolution of particle displacement Ax{t) during the recovery 

of a wormlike micellar solution after removing the trap, initially moving at different 
velocities v. From bottom to top: v = 1.60/J,ms“^ (+), 6.40p,ms“^ (*), 32.00^ms“^ 
(x). Time evolution of the normalized particle displacement [A — Ax{t)]/A for (b) 
a wormlike micellar solution, same data and symbols as in Fig. [H^a), (c) PAAm 
solution at u = 2.13/rms“^ (+), 4.00/rms“^ (+), 12.80/rms“^ (x), and (d) A—DNA 
at u = 1.30/rms“^ (+), 2.10^ms“^ (*), 3.20^ms“^ (x). The solid lines are hts to 
Eq. ID). 


on the time evolution of its x coordinate. We define ^2 = Os as the time at which the 
trap is turned off and Ax(t) = x(t) — x(t 2 ) as the distance travelled backward by the 
particle during the recovery of the fluid. For a hxed velocity v of the trap, the previous 
protocol is repeated at least 20 times in order to perform an ensemble average of Ax. 
We typically span values of v between 0.4/ims“^ and 40/ims“^. For these values of 
V the Deborah number is between 0.025 and 2.5, whereas the Weissenberg number is 
between 0.125 and 12.5, thus fully probing both linear and nonlinear viscoelasticity of 
the investigated fluids. 
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In Fig. we plot the time evolution of the displacement Ax{t) of a r = 1.6 fim particle 
after removal of the trap moving at different initial velocities n in a wormlike micellar 
solution (CPyCl/NaSal 5 mM, T = 303.16 K). We hnd that Ax{t) strongly depends 
on V, with a monotonic increase over a hnite time-scale until the complete relaxation 
of the system to thermal equilibrium, at which Ax{t) saturates to a constant value 
A = Ax{t —)■ cx)). In order to understand this dynamics, in Fig. Et^b) we represent in 
a semi-log plot the time evolution of the normalized displacement, [A — Ax{t)]/A. As 
a result all curves start at 1 and decay monotonically to 0 regardless of the value of v. 
We observe that all the curves exhibit an initial fast decay, followed by a slower one, 
which suggests the existence of two exponential relaxations. Indeed, we verify that for 
all V the time evolution of Ax{t) can be htted to the function 


Ax{t) = A — Afcxp 




( 1 ) 


with two relaxation times Tf and Tg and two amplitudes A{ and Ag which satisfy 
A = As + A{. This is shown in Figs. |2]^a) and EJ^b), where the hts to Eq. ([1]) are 
represented as solid lines. Furthermore, we hnd that this transient behavior of the 
particle position is not specihc to the recovery of wormlike micellar solutions. For 
example, in Figs. EJc) and EJd) we plot the time evolution of [A — Ax{t)]/A for a 
r = 1.6/im particle recoiling in a PAAm solution and in A-phage DNA, respectively, 
for different v. We hnd that, in these cases Eq. ([T]) also describes the time evolution of 
Ax{t), even when the rheological properties of such viscoelastic huids are more complex 
than those of wormlike micelles [33]. 

To support the existence of two distinct relaxation processes, we studied the 
dependence of the htting parameters of Eq. a on the velocity v. For instance, in 
Fig. El^a) we plot the two relaxation times Tf and Tg as a function of v during the 
recovery of a wormlike micellar solution (CPyCl/NaSal 5 niM, same sample as in 
Fig. Et^a)). We hnd that the values of these two time-scales are diherent by one order 
of magnitude, and interestingly, both are independent of v. For this reason, we will 
refer to them in the following as fast and slow relaxation times, respectively. Their 
corresponding amplitudes, Af and Ag, are plotted in Fig. El^b). We observe that, unlike 
the relaxation times, the amplitudes strongly depend on v. In particular, while Af 
increases monotonically with increasing v and then it levels oh, Ag has a maximum and 
then decreases with increasing v. 

For the values of v we span in our experiments, both linear and nonlinear 
microrheological response of the viscoelastic huid is induced by the particle motion 
previously driven by the tweezers. This is shown in Fig. |3]^c) where we plot the ehective 
viscosity p of the huid, measured by active microrheology during step (ii) of Fig. E^a), 
as a function of the imposed velocity v. As highlighted in the inset of Fig. |3]^c), at small 
f < 3 yums“^, 7] is constant and equal to the zero-shear viscosity rjo = 0.045 ± 0.005 Pa s 
(dashed line), because the strain rate imposed by the particle is sufficiently slow for 
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Figure 3. Dependence of (a) relaxation times Tg (□), T{ (o), and b) amplitudes Ag 
(□), A{ (o), obtained from the fit to Eq. (H)), on the driving velocity v. The symbols 
(*) correspond to A = As + Ai, whereas the dashed lines represent the values Tgg and 
Acfi computed from Eq. (0 for the rheological model ([2]). (c) Effective viscosity g of 
the fluid as a function of the driving velocity v determined by active microrheology 
during step (ii). Inset: log-log plot of 77 as a function of the corresponding Weissenberg 
number, same data as in main Figure. The dashed line represents the value of the zero 
shear viscosity 770 . 


the fluid to lose its memory on previous microstructural deformations. Consequently, 
Newtonian-like behavior is observed through the validity of the linear-response Stokes 
law: fd = QurrjQV. On the other hand, at sufficiently large v > 3pms“^, the particle 
motion is able to induce an orientational order of the fluid microstructure. As a 
result, the fluid exhibits thinning, where p decreases dramatically with increasing v en- 
Note that the emergence of nonlinear rheological effects of the fluid due to the local 
deformation created by the particle are better quantified by means of the Weissenberg 
number Wi. In our system, Wi can be estimated by the product of the characteristic 
rate of deformation ^ ^ times the largest relaxation time of the system Tg; 

Wi = VTs/{2r). Then, nonlinear non-Newtonian behavior is expected at Wi > 1, at 
which the local microstructural deformation of the fluid is so strong that the Stokes 
law must break down. Indeed, in Fig. El^c) we show that for the spanned values of 
V, the effective viscosity of the fluid exhibits the aforementioned transition at Wi ~ 1 
and therefore we cover both linear and nonlinear microrheological behavior. Hence, our 
results demonstrate that the double-exponential relaxation of the particle position is 
a rather general process during the recovery of viscoelastic fluids with a well-defined 
zero-shear viscosity, regardless of the value of the Weissenberg number. 

We point out that the dependence on the velocity v of the relaxation times Tf 
and Tg, and the amplitudes Af and Ag is not particular to the particle motion during 
the recovery of wormlike micelles, but also found for the other investigated viscoelastic 
fluids. In Figs. Il]^a)-(c) we plot the results for an aqueous semidilute polymer solution of 
PAAm (Mu, = 18 X 10® at 0.05 % wt, temperature T = 293.16 K), whereas the results 
for A-phage DNA (500 /ig/niL in 10 mM Tris-HCl and 1 niM EDTA, temperature 
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Figure 4. Dependence on imposed velocity v of (a) relaxation times Tg (□), T{ (o), 
(b) amplitudes Ag (□), Af (o), A = As + Af {*) and (c) effective viscosity y (o) for a 
r = 1.6 ^m particle in a polymer solution of PAAm. Dependence on v of (d) relaxation 
times Tg {□), Tf (o), (e) amplitudes Ag (□), Af (o), and A = Ag+Af (*) and (f) effective 
viscosity y (o) for a r = 1.6 ym particle in A—DNA. The dashed lines in Figs.lU^a) and 
md) correspond to Teff whereas in Figs.|Ub) andlUJe) they represent Aeg. See text for 
explanation. 


T = 298.16 K) are displayed in Fig. |4]^d)-(f). Once again, we observe that the two 
distinct relaxation time-scales remain independent of v both in the linear and nonlinear 
response regime, whereas the amplitudes exhibit a very similar dependence to that 
observed for the wormlike micellar solution. Thus, our findings suggest that only the 
amplitudes of Ax{t) and not its relaxation times encode rheological information of the 
local deformation of the fluid microstructure created by the particle, the latter defining 
the time-scales at which the system retains memory on the energy stored by such an 
elastic deformation. 

4. Discussion 

In order to gain an understading of this obviously generic behavior of Ax{t), we now 
investigate whether a Langevin model can reproduce the particle dynamics after sudden 
removal of the driving force after dragging the particle over a large distance Xq > 2r 
















Transient dynamics of a colloidal particle driven through a viscoelastic fluid 


9 





Figure 5. (a) Time evolution of the coordinates x (red solid line) and y (black 

dashed line) of a r = 1.6/rm particle in a wormlike micellar solution in response to 
a small-amplitude sinusoidal motion of the optical trap (fc = 2 x 10“®Nm“^) with 
amplitude Xq = 200 nm and angular frequency = 27r rad s“^. (b) Experimental 
power spectral density of x (dotted lines) at different oscillation frequencies fl/(27r) 
of the trap. The peaks at each imposed Q,/{2'k) are indicated by circles. The dashed 
line represents Eq. whereas the solid line corresponds to the total power spectral 
density described by Eq. ®. (c) Schematic representation of the particle dynamics 
dragged at constant velocity —v through a fluid with linear shear modulus Q (upper 
panel) and during the relaxation after removing the driving force (lower panel). See 
text for explanation. 


through the fluid, i.e. after inducing a large strain Xo/(2r). For this purpose we 
consider the simplest case of a particle driven by the moving trap in a wormlike micellar 
solution. In general, under sufficiently small applied shear strain rate 7 , the resulting 
shear stress of a viscoelastic fluid a is determined by the linear shear modulus G\ 
a{t) = G{t — t')'j{t') dt'. In the case of a wormlike micellar solution, G can be 
simply characterized by a single stress relaxation time r, the viscosity of the solvent 700 
and the zero-shear viscosity rjo (Jeffreys model) [M] 

G{t) = 2T]^6{t) + exp • (2) 

Therefore, according to commonly accepted assumptions, the steady-state dynamics 
of the X coordinate of the particle is described in the linear-response regime by the 
non-Markovian Langevin equation [53] 

f r{t — t')x{t') dt' = —k[x{t) — Xo{t)] + ^{t). 


(3) 
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In Eq. ([3]), r{t — t') = GirrGfl — t') is a memory kernel related to the shear modulus G of 
the fluid given in Eq. (|2]), whereas is a coloured Gaussian noise with autocorrelation 
function {^{t)^{t')) = fc^TEdt — t'\). 


4 . 1 . Transient motion under small deformation 


We hrst experimentally verify that the Langevin Eq. (|3]) correctly describes the 
equilibrium dynamics of the particle and its linear response to small-amplitude 
oscillatory forces exerted by the optical trap. This implies that, for sufficiently small 
deformations of the fluid, bulk rheological properties can be unambiguously determined 
by measuring the position of the embedded particle, as widely assumed in passive and 
linear active microrheology. For this purpose, we apply to the particle a sinusoidal 
force kxoit) by moving the trap according to Xoifl) = Xosin(r2t), with amplitude 
Xq = 200 nm 2r at different driving frequencies Q/{27r). The characteristic strain 
of the fluid locally induced by this oscillatory motion is Xo/{2r) = 0.06. The typical 
time evolution of the particle position in response to such a driving force is illustrated 
in Fig. ISt^a). We observe that, while the coordinate x responds to the force at the same 
excitation frequency r2/(27r), the coordinate y remains unaffected. Therefore, linear 
response dynamics can be tested by means of x, whereas thermal equilibrium properties 
can be analysed by means of y. This can be done, e.g. by comparing the experimental 
power spectral densities of such stochastic signals to those predicted by the Langevin 
model (l3j). For a sinusoidal driving force, Eq. ([3]) yields the following expressions for 
the power spectral densities of x and y in terms of the linear storage G'(/) and loss 
G"{f) modulus of the fluid 


{\x{f)\^)=S,asif) + SaM), ( 4 ) 

mm = S,as{f), (5) 


where 

Q (n = ‘^ML _ 67rrG"(/) _ 

71 f [/c-|-67rrG'(/)]^-|-[67rrG"(/)]^ ’ 

. (r. Kf-B 

^ 2 [k + 67rrG'(/)]2 + [67rrG"(/)]2 ’ 


( 6 ) 

( 7 ) 


are the contributions to the power spectral density due to thermal fluctuations and to 
the active driving force, respectively. In particular, the prefactor ksT in Eq. ([6]) results 
from the thermal origin of the fluctuations, while the delta function d{f — ^) in Eq. ()7|) 
mirrors the mechanical response of the system at frequency / = Indeed, we hnd 
the expressions dl])-© derived from Eqs. ([2]) and ([3]) correctly describe the measured 
power spectral densities of x and y for all the investigated frequencies. This in turn 
validates such a Langevin approach in the linear response regime, provided that the 
particle performs small displacements in the fluid, i.e. when sufficiently small flow and 
deformations are induced. We show this in Fig. |5]^b), where we demonstrate that, at the 
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unexcited frequencies / 7 ^ the power spectral density is correctly described Eq. (j6|) 
(dashed line), where G\f) and G”{f) are obtained from the Fourier transform of Eq. ([2]). 
On the other hand, an excellent agreement between the experimental peaks of (|£(/)P) 
at the excited frequencies / = ^ (circles) and Eq. (|1]) (dashed line) is observed. 

Once we have verihed that the Langevin model ([3]) correctly describes the particle 
dynamics under sufficiently small deformations of the fluid, we now focus on its transient 
behavior after removal of the driving force. A schematic representation of the average 
forces acting on the particle is shown in Fig. |5](c). Interestingly, Eq. (|3]) predicts that 
after dragging the particle at velocity —v at time t < 0 through a medium with linear 
viscoelasticity (EJ, it recoils with an initial velocity {rjo — rjoo)v/r}oo at time f = 0 upon 
removal of the trap. Moreover, the subsequent displacement Ax{t) evolves in time f > 0 
according to a single exponential relaxation 


Ax{t) = Aeff 


1 — exp 


with amplitude and relaxation time 



( 8 ) 


Aeflf — 


1 _ Ha. 

ho 


nr, 


hoc 


(9) 


respectively. See [Appendix A| for a more detailed derivation of Eqs. ([8]) and ([9]). We 
notice that, according to Eq. ([2]), the amplitude Agff depends linearly on v and encodes 
information of the rheological properties of the fluid. Indeed, the driven motion of the 
particle induces a strain in the fluid, which is captured in the linear rheological model 
of Eq. (E]) as an effective elastic compression ~ vr, as depicted in the upper panel of 
Fig. |5](c). This results in a total displacement of the particle proportional to v upon 
removing the driving force, as schematized in the lower panel of Fig. El(c). On the 
other hand, the relaxation time Tes of the particle displacement is independent of v and 
always smaller than the stress-relaxation time of the fluid r. This is due to presence 
of the solvent, whose viscosity r^oo makes the particle become at rest in the effective 
Langevin model ([3]) before the fluid completely relaxes its stress initially stored by the 
elastic compression induced during the active driving process. 


4-2. Transient motion under large deformation 

We now compare the values of the effective parameters Tes and Aes with those we obtain 
experimentally for the colloidal particle after performing large displacements Xq > 2 r 
through the fluid, i.e. after inducing a characteristic strain Xo/(2r) = 5. In Fig. [3](a) we 
show as a dashed line the value of the effective relaxation time Teg, where the parameters 
ho 5 hoo and r are determined by passive microrheology during step (i) of Fig. [I](a). We 
observe that the values of the fast relaxation time Tf obtained from the ht to Eq. o 
(circles) agree well with Tes, which implies that the fast relaxation process mirrors the 
effective elastic relaxation of the probe damped by the solvent viscosity. Nevertheless, 
in Fig. ini(b) we show that Aeg > Af for all n > 0 which suggests that the fast relaxation 
is just one of the possible channels for the complete relaxation of the system to thermal 
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equilibrium. As a matter of fact, in Fig. [3]|^b) we show that at small v, A = Ag + Af agrees 
well with Aeff. Thus, A actually quantihes the elastic compression of the fluid induced 
by the driven particle in the linear response regime. Systematic deviations of A from Aes 
show up with increasing v, though, because nonlinear microrheological response of the 
fluid is induced during the previous driven motion of the particle. In such a case, Eq. o 
fails to describe the viscoelastic properties of the fluid. We conhrm this in Fig. [3|^c), 
where we show that at small v ^3 /ims“^, i.e. Wi < 1, linear microrheological response 
is induced, i.e. t] = rjo, whereas at sufficiently large v > 3 /ims“^, i.e. Wi > 1, the fluid 
microstructure responds nonlinearly and the system exhibits thinning. Note that the 
values of v for which thinning emerges agree well with those for which large deviations 
of A from Agg are observed. Nonlinear microrheological behavior of the fluid affects 
separately the behavior of Af and As as well. In particular, Af exhibits a saturation 
whereas Ag has a maximum and then decreases for the values of v at which prior thinning 
is observed. These results clearly demonstrate that nonlinear response of a viscoelastic 
fluid, induced by active microrheology, gives rise to nonlinear transient behavior of the 
probe upon removal of the active driving. 

We also check that the previous mechanism applies to the recovery of the polymer 
solution and of A—DNA. In absence of simple rheological models for the linear shear 
modulus G(t) of such viscoelastic fluids, in analogy to the observations for the wormlike 
micellar solutions, we hrst assume that Tes = Tf and t = Tg, which implies that 
Voo = Vo'Tes/T = rjoTf/Tg. Then, accordingly to Eq. ([2]) we compute the effective 
amplitude by means of 

Agff = u(rs - Tf). (10) 

Surprisingly, we observe in Figs. 111(b) and 111(e) that the values of Ag^ obtained from 
Eq. firoj) (dashed lines) agree very well with the total amplitude A = Ag + Af (asterisks) 
at sufficiently small velocities v for both fluids, thus making evident the robustness of the 
two distinct time-scales T{ and Tg regardless of the details of the rheological properties 
of the embedding viscoelastic fluid. On the other hand, for velocities at which thinning 
is induced, i.e. at which rj decreases with increasing particle velocity, the amplitudes 
exhibit a dependence very similar to that observed for wormlike micelles. While Af 
saturates with increasing v, As has a maximum and then decreases. The consequence 
is that the total amplitude A levels off with increasing velocity, thus signalling the 
emergence of a nonlinear regime where the elastic compression of the fluid induced by 
the particle becomes independent of v. Hence, our results demonstrate the generality 
of the transient microrhelogical behavior for very distinct viscoelastic fluids both in the 
linear and nonlinear response regimes. 

The previous Endings make evident that the transient rheological behavior of 
viscoelastic fluids probed by microrheology does not only depend on the characteristic 
rate of local deformation v/{2r), but also on the total strain Wo/(2r). On one hand, 
for small Ao/(2r) 1, the Langevin model ([3]) correctly describes the dynamics of the 

particle either in thermal equilibrium with the surrounding fluid or externally driven into 
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t [s] x[s] 

Figure 6. (a) Example of time evolution of imposed stress and (b) resulting strain 
during a shear-flow cessation experiment of a micellar solution (CPyCl/NaSal 7 mM) 
in a Couette geometry. The dashed line in Fig. IH^b) represents an exponential fit with 
relaxation time r = 1.8 s. (c) Comparison between the slow relaxation time Ts of 
various viscoelastic fluids determined by the transient motion of particles of different 
radii r and the equilibrium stress-relaxation time r. The symbols are: CPyCl/NaSal 
5 mM, r = 0.9/rm (x), r = 1.3/rm (o), r = 1.6/iin (*), r = 2.2/im (•); A—phage 
DNA, r = 1.6/xm (x); PAAm, r = 1.6/xm (□); CPyCl/NaSal 7 mM, r = 1.3/xm (<), 
r = 2.2 yvci (i>); CPyCl/NaSal 9 mM, r = 1.3/xm (v), r = 1.6/xm (a). Inset: expanded 
view of main figure around 0.9 s. 


linear response regime by small-amplitude forces, where vanishingly small deformations 
of the fluid are induced. On the other hand, a different transient behavior is observed for 
sufficiently large strain Xo/(2r) > 1, where such a Langevin model ceases to be valid. In 
particular, it fails to reproduce the slow relaxation inferred from Eq. ([2]) after cessation 
of a large deformation imposed by the driven particle. The reason is that Eq. (E]) only 
accounts for the coarse-grained drag force exerted by the fluid. Such a viscoelastic drag 
force is enough to describe equilibrium thermal fluctuations of the particle position 
without recourse to any information on the surrounding flow held. However, it does not 
incorporate the effect of the transient how and deformation localized around the particle 
when dragged over large distances, and which cannot be neglected during the recovery 
process of the huid. In the simplistic Langevin picture, the only possible channel for 
the system to reach thermal equilibrium is by means of the ehective relaxation of the 
particle position. While this relaxation probes the linear bulk properties of the huid, it 
ignores the complete stress relaxation of the local viscoelastic how held. 

In fact, the aforementioned drawbacks suggest that, in reality, the residual slow 
relaxation represents a second mechanism to completely release the stress stored by 
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the strained fluid. The two relaxation mechanisms allow in turn the fluid to reach 
thermal equilibrium in conjunction with the particle motion. These observations hint 
at the stress-relaxation time r of the fluid as the slow relaxation time Tg. In order to 
test this hypothesis, we perform using a rotational rheometer (HAAKE RheoStress 1, 
double-gap cylinder geometry) the macroscopic version of the transient microrheology 
experiments, This is illustrated in Figs. | 6 ])a) and[ 6 ])b), where we plot an example of the 
time evolution of the shear stress a applied to a wormlike micellar solution, and the 
resulting shear strain 7 , respectively. After reaching a steady state with constant shear 
rate 7 at time t < 0, the shear stress is suddently removed at f = 0 , and afterwards 
7 decays in time, until the fluid completely recovers. Then, we extract the relaxation 
time r of the fluid from the exponential fit 7 (f) = 7 oexp(—f/r). In Fig. [H](c) we plot 
Tg, as a function of the corresponding r for all the investigated viscoelastic fluids. We 
And that within the experimental errors, both time-scales are equal. In Fig. | 6 ])c) we 
also verify that the relaxation times Tg are independent of the size of the probe particle. 
See Appendix B for more detailed results obtained for different particle sizes. The 
independence of Tg from v, r and Xq, i.e. from the detailed flow and strain field around 
the particle, demonstrates that it is a physical property inherent to the viscoelastic fluid. 
Thus, this time-scale Axes a value above which the memory on the storage of elastic 
energy by the fluid microstructure is lost. Hence, the the slow transient dynamics of a 
colloidal probe after removal of a driving force represents a new straightforward method 
to measure the stress-relaxation time of viscoelastic fluids, provided that the previous 
driven motion induces a sufficiently large strain. This can be directly carried out by 
means of particle tracking without the need of any model for the dynamics or any further 
elaborate Fourier or Laplace analysis, as commonly done in passive and small-amplitude 
active microrheology [371138] . 

We point out that the slow relaxation is most prominent close to the onset of 
nonlinear microrheological response around Wi ~ 1, as verified by the presence of the 
maxima of Ag in Figs. [21(b), 111(b) and 111(e). Therefore, our findings have important 
consequences on the interpretation of transient behavior in active microrheology under 
large deformations [201 ISHl ES]. This is because they provide experimental evidence 
that, even at small Weissenberg number, where only linear microrheological response 
is expected to occur, the response of the colloidal particle is strongly affected by the 
relaxation of the local flow and strain field around it. This is in contrast to passive and 
active microrheology under small-amplitude deformations, where only bulk properties of 
the fluid are probed by the particle. In addition, our results have significant implications 
for the behavior of micron-sized objects in transient viscoelastic flows, such as those 
typically found in microchannels, where large strains are commonly induced by high 
flow rates and sharp constrictions. In such systems, unlike the case of Newtonian 
liquids, the conventional Langevin description for the suspended particles involving only 
coarse-grained drag forces and noise must break down. Furthermore, our results are 
particularly important for the understanding and modelling of the motion of biological 
microswimmers, e.g. bacteria and spermatozoa, whose natural environment is frequently 
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viscoelastic ilElIID]. Such active systems usually operate at high Weissenberg number 
and travel large distances through the viscoelastic medium, thus inducing complex 
flow flelds and large strains. Therefore, the relaxation of the elastic energy stored by 
the surrounding fluid microstructure can strongly impact their swimming mechanisms 
through the different processes found in our experiments. 

5. Summary and conclusion 

We have experimentally studied the motion of a colloidal probe during the local strain 
recovery of several viscoelastic fluids upon shutoff of an active driving force. We have 
shown that, even when a Langevin model correctly describes linear response dynamics 
of the particle position under vanishingly small strain, it cannot capture its transient 
dynamics during the recovery process when large deformations are previously induced 
by large displacements of the particle through the fluid. In particular, we have found 
that this transient dynamics proceeds via a double exponential decay, which reveals 
the cooperative relaxation to thermal equilibrium of the particle and the surrounding 
strained fluid. We have shown that the fastest relaxation mirrors the viscous damping 
of the particle by the solvent, whereas the slow one results from the relaxation of the 
viscoelastic matrix. Thus, the two relaxation processes represent dissipation channels 
for the elastic energy stored by the fluid microstructure upon cessation of a sufficiently 
large deformation. We have demonstrated that, while the amplitudes of these relaxation 
modes encode either linear or nonlinear microrheological information of the fluid induced 
during the prior driving, their relaxation times are insensitive to it. Consequently, this 
transient microrheological method allows to determine unambiguously stress-relaxation 
times of micro-litre samples of viscoelastic fluids. 
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Appendix A. Transient particle’s dynamics in a Jeffreys fluid 

For the Langevin model of a particle moving in a Jeffreys fluid, i.e. Eqs. ([2]) and ([2]), an 
analytical expression for the time evolution of its position Ax{t) after cessation of an 
externally applied force can be readily derived. This can be achieved by realizing that, 
for r]o > r]oo > 0 and r > 0, the non-Markovian equation of motion ([3]) for x can be 
written as a linear system of two Markovian Langevin equations [31] 


x{t) 




k[x{t) - Xo{t)] ^ / ksT 


6vr?7oor V ?,'Kri^r 
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u{t) = -[x{t) - u{t)] + 

T 


knT 


-ut). 


(A.l) 


3vr(r7o - y]oo)r 

where Ci and C 2 are white noises of mean {Q) = 0 and correlation {Ci{t)Cj{t')) = 
Sij6{t — f), i,j = 1, 2, whereas the auxiliary variable u is dehned as 

1 rt 


u{t) = - 

T 


exp 


t-t' 


x{t') + Ti 


ksT 


3tt{t]o - rioo) 


Ut') 


dt'. (A.2) 


Without loss of generality, we dehne the particle position at time t = 0 at which 
the external force is removed {k = 0) as a;(0) = 0. Therefore, the solution of the linear 
system flA.ip at subsequent times f > 0, averaged over a large number of realizations of 
the noises and ( 2 , can be expressed in terms if the initial conditions x(0), m( 0) and 
'u(O) at t = 0 as 


{x{t)) 

{u{t)) 


■ [r\\hoo 

xfO)—r 
Vo 


■u(O) + hfO)—r 
Vo 


1 — exp 



L V 

VO / J 


1 — exp 



L V 

VO / J 


(A.3) 


In order to obtain the initial conditions of Eq. flA.3p . we take into account that 
at f < 0, the particle is dragged in a steady state by the optical trap at velocity —v, 
i.e. {x{t)) = —vt. On the other hand, upon removal of the external force at f = 0, the 
ensemble average of the viscoelastic drag force Fdrag{t) = T[t — t')x{t') dt' in Eq. ([3]) 
is zero. Therefore, after integration by parts and using the auxiliary variable u dehned 
in Eq. (1A.2D . we hnd the relation 


0 = {Fdrag{d)), 

= Gvrr 


• /'n^ ho Voo 

VooX^)) -m(0) 


r 


(A.4) 


where 


' 1 .( 0 ) = 


r j 

= VT. 


(a:(f'))exp ( - ) dt', 


T 


(A.5) 

Eqs. flA.4p and flA.Sp yield the following expressions for the initial conditions of 
Eq. flA.3p in terms of the initial driving velocity v and the rheological quantities of the 
Jeffreys model 


/p|\ PI • /p|\ Vo Voo 
xfO) = 0, x(0) =- V, 

Voo 

m(0) = vr,u{0) = —V, 


(A,6) 


Finally, by inserting flA.6p into Eqs. flA.3p . we hnd the analytical expression of Eq. ([3]) 
for the distance recoiled by the particle Ax(t) = {x(t) — x(0)) after removal of the 
imposed velocity. 
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Figure Bl. Dependence on imposed velocity v of (a) relaxation times Tg (□), Tf 
(o), (b) amplitudes Ag (□), A{ (o), A = Ag + Af (*) and (c) effective viscosity g (o) 
for a r = 0.9 ym particle in a wormlike micellar solution. Dependence on v of (d) 
relaxation times Tg (□), Tf (o), (e) amplitudes Ag (□), Af (o), and A = Ag + Af {*) 
and (f) effective viscosity rj (o)for a r = 1.3/rm particle in the same micellar solution. 
The dashed lines in Figs. IBlf ab IBlf b'). IBlf d! and IBlf ei correspond to Tgff and Aeff 
described in the main text, whereas those in Figs. IBlT ci and Figs. IBlT fi represent the 
zero-shear viscosity go of the fluid. 


Appendix B. Particle size 

In this section, we show the resnlts for the transient dynamics of particles of 
different radii r dnring the microstrnctnral recovery of a wormlike micellar solntion 
(CPyCl/NaSal 5 mM at T = 303.16 K). Figs. IBlf al-fcl are the results for r = 0.9/im 
whereas those for r = 1.3 pm are plotted in Figs. IBir dl-lfh We point that the 
dependencies of the two relaxation times Tf and Tg, (Figs. IBlf al and IBlf dl). the 
amplitudes Af, Ag and A = Af + Ag (Figs. IBlf bl and IBlf eA and the viscosities rj 
(Fig. IBlf cl and IBlf fA are very similar to those shown in Fig. [3] the main text for the 
same micellar solution using a particle of radius r = 1.6 pm. In particular, the values of 
Tf and Tg do not depend on the particle size r, whereas the amplitudes are independent 
of r in the linear response regime at which p = po- The only difference is observed in the 
nonlinear values of the amplitudes for which thinning is previously induced at sufficiently 
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large v. In such a case, Figs. IBlT bi and IBlIf ei suggest that with decreasing particle 
radius, the total amplitude A decreases. For instance, A 1.0 /rm for r = 0.9 /im, 
A 1.5 fim for r = 1.3 /im, and A ^ 1.7 fim for r = 1.6 /im. This can be explained by 
the fact that, at hxed v and Xq, the smaller the particle radius r, the stronger the local 
deformation rate v/{2r) and the total strain Xo/(2r) of the fluid induced by the active 
driving, thereby resulting in a weaker nonlinear strain recovery. 
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